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Abstract 


The finite element presented in this report is useful for analyzing plates and 
shells by the direct stiffness method. Both conventional thin plates and shells, 
as well as sandwich structures, are included. The stiffness matrix derivation is 
both simple and straightforward, since it is based on previously published ma- 
terial. A few examples are included to illustrate the accuracy obtained with this 
element. 
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Stiffness Matrix for a Triangular Sandwich 
Element in Bending 


I. Introduction 

This report presents various phases of the development 
of a triangular finite element in bending. A sandwich 
element is chosen because it provides the simplest basis 
for obtaining a satisfactory stiffness matrix. At the same 
time, the sandwich element is applicable to sandwich 
panels, as well as to conventional thin plates and shells. 

The triangular element stiffness matrix, plus the nu- 
merical proceedings of the stiffness method, yields de- 
flections and stresses for the most general classes of plate 
and shell problems. Although results will be approximate, 
experience has shown that they are usually adequate for 
engineering design purposes. 

Interest in the triangular element in bending during 
the late 1950’s, led to work directed toward finding a 
suitable stiffness matrix. These efforts proceeded along 
two totally different paths. The first considered the ele- 
ment as a typical thin-plate unit. By choosing a deflection 
function, w ( x , y), straightforward concepts of the stiffness 
method could be followed in deriving the element stiff- 
ness matrix. References 1-4 give examples of this early 
work. Results were generally unsatisfactory and could 
not be used in engineering design. The difficulty resided 
in the selection of a satisfactory deflection function. 


At about the same time, R. J. Melosh, then at The Boe- 
ing Company, recognized that a second path was avail- 
able. The stiffness matrix for the triangle in plane stress 
was already known (Ref. 5). Hence, it was decided to 
use these parts as the cover sheets of a sandwich element 
and, on that basis, obtain a stiffness matrix for the tri- 
angle in bending. Melosh (Ref. 6) in his paper states that 
“publication was delayed due to dissatisfaction with the 
derivation in its initial form.” Nevertheless, a correct result 
for the sandwich element had been found. This led to 
satisfactory data when applied to plate and shell prob- 
lems. Reference 7 gives a number of examples, the solu- 
tions to some of which were reached by using the original 
stiffness matrix obtained for the triangular sandwich 
element. 

In the meantime, it has been recognized that the tri- 
angular element in bending presents difficulties which 
had not been encountered in the initial work on finite 
elements. By adopting the standard procedure of simply 
considering nodes to exist at the three vertices of the 
triangle, it was eventually realized that a satisfactory func- 
tion w ( x , y) could not be established. Fraeijs de Veubeke 
first suggested the need for introducing additional nodes 
in such cases, and this led to the so-called HCT (Hsieh, 
Clough, Tocher) triangle as presented by Clough and 
Tocher in Ref. 8. By introducing an additional node at 
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the centroid of the triangle, it became possible to replace 
the original area with three subtriangles. This, in turn, 
permitted much greater generality to be exercised in 
forming w (x, y). Clough and Tocher show in Ref. 8 that, 
on this basis, a greatly improved bending stiffness matrix 
for the conventional triangle in bending could be derived. 

The sandwich element discussed in this report leads to 
the same stiffness matrix obtainable from Ref. 6. The 
same element is used in the COSMOS program at The 
Boeing Company and in the SAMIS program developed 
at the Jet Propulsion Laboratory, However, the derivation 
is entirely different from that given in Ref. 6. The devel- 
opment of the stiffness matrix proceeds in a straightfor- 
ward manner if based on a clear statement as to the nature 
of the element. An explicit final form is obtained for the 
stiffness matrix. This, in turn, proves to be useful for 
looking into certain questions arising with respect to this 
element. 

It # is important to recognize that the sandwich-type 
element enjoys some unique characteristics which make 
it useful to the structural designer. Among these are: 

1. Orthotropic, or even aeolo tropic, material properties 
may be represented. Furthermore, these properties 
need not be the same for cover sheets and connect- 
ing shear core. 

2. Design orthotropy, as illustrated by corrugation or 
stiffeners held between cover sheets, may also be 
represented by this element. 

3. Plate or shell sections having appreciable shear flex- 
ibility are reasonably idealized by the sandwich 
element. 

4. Structures containing plates or shells, reinforced 
with beams or frames, present idealization problems. 
In general, finite plate and beam elements will not 
remain compatible along common edges as deflec- 


tion takes place. Since the sandwich plate and beam 
elements do not undergo curvature, the opportunity 
for achieving such compatibility is greatly enhanced 
if the sandwich elements are used. 

5. The triangular sandwich element preserves both 
displacement and ‘slope” continuity across edges 
common to adjacent triangles. Furthermore, this is 
true even when the triangles lie in different planes, 
as when used to represent a curved shell. 

6. Since this element remains flat when representing 
plate or shell bending, it may be used with the initial 
stress stiffness matrix developed by the author 
(Ref. 9) to represent large deflection and stability 
behavior. 

II. The Isotopic Triangle in Bending 

A. Sandwich Element— Description 

The triangular element presented herein is visualized 
as consisting of identical top and bottom cover sheets, 
separated by thin shear webs located at each edge of the 
triangle (Fig. 1). 

To specify the element, it is necessary to give its over- 
all thickness d, cover sheet thickness h, and shear web 
thickness t w , as well as E, G, and v , in terms of the 
actual plate. The basis for doing this will be taken up 
prior to deriving the element stiffness matrix. 

Bending moments are resisted solely by the cover 
sheets. These plane stress members are restricted to a 
constant state of strain. As such, their stiffness matrix is 
given by Eq.(B-3), Ref. 5. Shear webs transmit the trans- 
verse loading. Each is assumed to be in a state of pure 
shear. As such, their stiffness matrix is given by Eq. 
(4.85), Ref. 10. The stiffness matrix for the total element 
will be developed by using these component stiffness 
matrices. 



CONNECTING SHEAR WEBS LOCATED 
ALONG EACH EDGE OF TRIANGLE 

TOP COVER SHEET (THICKNESS=A ) 

BOTTOM COVER SHEET (THICKNESS * A) 


Fig. 1. Component parts of the triangular sandwich element 
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The corresponding sandwich beam element is discussed 
in Ref. 10. Reference 10 also shows how well such an ele- 
ment will predict actual beam behavior. Both the beam 
and triangular elements remain flat as they deform under 
loading. Nevertheless, the ability of the sandwich element 
to predict beam behavior provides encouragement for 
extending this concept to the more complex case of the 
i triangle. 

It should be kept in mind that the top and bottom 
covers are strained asymmetrically. For the beam element 
! this means that, if the upper flange is in compression, the 

lower flange will experience an equal tensile strain. 

B. Sandwich Element-Properties 

I. Displacements and forces. If it is assumed that the 
element lies in the xy plane, then displacements u and v 
at each node define the deformation in the cover sheets. 
Normal displacements tv define the shear web deforma- 
tion at each node. It will be convenient to replace u and 
v by slopes d x and 6 y for the complete element. Conse- 
quently, three nodal displacements exist at each node, 
leading to a 9 X 9 stiffness matrix for the element. A 
sketch of the element is shown in Fig. 2a. The middle 
plane is shown in Fig. 2b. 



Fig. 2. Nodal representation for the sandwich element 


Asymmetrical behavior of the top and bottom cover 
sheets are defined by 

Ua = — Ua’ Va ~ — Va' (la) 

and 

Wa ~ Wa’ = tVi (lb) 

where 

a = a, b,c 


and 


i = 1,2,3 


The slopes may be introduced in place of u and v by 
referring to Fig. 3 which demonstrates that the following 
equations apply: 


Va 

d/2 9 


** = 


ti a 

1/2 


( 2 ) 




(b) 


Fig. 3. Relation between cover sheet displacements 
(u, v) and corresponding sandwich element 
displacements i0 x ,0 v ) 


The signs in Eq. (2) follow from Fig. 2b and Figs. 3a 
and 3b. 

Writing Eq. (2) at all nodes and collecting results 
gives 



(3a) 


or simply 

u = ^T0 (3b) 

Note: Boldface symbols represent matrices. 

All elements not shown in T, Eqs. (3a) and (3b), are zero. 

Cover sheets carry in-plane stresses a x , <r y , and T*y. In 
deriving the cover sheet stiffness matrix these stresses 
are replaced by equivalent nodal forces, X a and Y«. Fig- 
ure 4 illustrates these nodal forces. 
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b 


(a) (b) 

Fig. 4. Nodal forces for cover sheets 


2. Cover sheet thickness . The cover sheet thickness is 
established by requiring a unit width of the sandwich 
panel to have the same moment of inertia as a unit width 
of the true plate. Or, from Fig. 5 

I (sandwich) = I (actual structure) 



The following relations for the forces are in harmony 
with Eqs. (la) and (lb) for displacements. 


In the last equation the thickness, d, of the sandwich may 
be selected arbitrarily. For convenience, put 



d — t 


(7a) 


(4a) Then 


and 



(7b) 


Z a = Z a > - Zj/2 (4b) 


The moments applicable to the sandwich element are 
then 


M Xl = — Y u d , M Vl ^X n d (5) 


Again, writing the last equations for each node and col- 
lecting results 


Equations (7a) and (7b) define the sandwich element cover 
sheets in terms of the actual structure plate thickness. 




Fig. 5. Unit strips of sandwich and solid plates 


3. Shear web thicknesses . The depth of the shear web 
is t, as given by Eq. (7a). To determine the appropriate 
thickness for the shear webs, both the solid triangle and 
the corresponding sandwich element will be subjected to 
pure shearing displacement modes. The sandwich ele- 
ment will then be required to develop the same energy 
as the solid triangle. One such mode is represented by 
w 1 0, while all other nodal displacements are zero. 
This mode is shown in Fig. 6. 


where elements not shown in the square matrix of Eq. 
(6a) are again to be taken as zero. Comparing Eq. (6a) 
with Eq. (3) leads to 

M — d T T X (6b) 

where T T is the transpose of T. 


3 
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If Fig. 6 is initially assumed to represent shearing of the 
true, solid element, the shear strain y 1 may be written as 



Hence 


yi 



U> i 


where A w is the shear web cross-sectional area. Repeat- 
ing this calculation for the other two modes and super- 
imposing gives the total energy C7 n as 


Un = 




( 9 ) 


where 


A = P a L 2 - 3 / 2 = area of triangle 


Strain energies U i and U u can now be made equal to 
each other by requiring the coefficients of the squares of 
the displacements to be equal. Doing this and solving for 
shear web areas gives 


Strain energy l/ 7 (1> associated with this deformation may 
then be written as 

u t v = jf J jrytdxdydz 

where 



t — G yi 

Hence 

v<» = — r ^lu>l 

Repeating for the two additional shearing modes and 
superimposing gives the total strain energy U { for the 
solid plate as 

( L U w l + L U w \ + L Li W V ( 8 ) 

Next, assume Fig. 6 to represent the sandwich element. 
Shear webs located along edges 1-2 and 1-3 are strained 
in this mode. In this case, the respective shear strains are 

u>i _ to 2 

Yi ~ 7 — Y2 “ 7 

^ 1-2 *-* 1-3 

Corresponding strain energy then takes the form 



Equation (10) specifies the shear web cross-sectional 
areas. The sandwich element is now completely defined 
such that it will represent the actual thin plate which it 
is intended to represent. 

C. Stiffness Due to Cover Sheets 

The cover sheet stiffness matrix of Eq. (B-3), Ref. 5, 
relates X«, Y« nodal forces to the corresponding v a 
nodal displacements. For the top cover, the stiffness 
equation may be expressed as 



or 

X - K cs u (lib) 

where, from Ref. 5 
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Ua V a U}> V b U c V c 



A!X 2 3 + y\ 3 


(SYMMETRiC) 


^ 2 ^ 32^23 

*23 + j 


Eh<t> 

A.l* 23*31 + 

A,x 13 t /23 + >^ 23^13 1 A x x 2 3 + y\ 3 


2 

A1X231/13 a vx 13 y 23 

\ 

* 23*31 d" Aiy 2 3*31 [ A 2 Xi3t/ 3 i 

+ ^y 2 3 ; 


a 1X1 2X23 + y 12^/23 

1 

A-lX 2 lt /23 d" * / * 32 {/l 2 \ AiXi 2 *31 “f~ t/l 2*/31 

A,x 12 t/,3 + 1*13^12; A!*j 2 + y'i-.\ 


_AiX 32 y,2 + vx 2 iX 2 3 

1 

*12*23 d" Ait/i 2 t /23 | AiX 13 t / 12 d~ VXi 2 tf 13 

^ 12^31 A A,J/, 2{/31 I A2X21J/12 1*12 ^ 1 ^/l 2 _ 


in which (11°) 


X «0 -Xa — Xfl, etc. 


1 

* 2(1 - v 2 ) A 


The forces and displacements of Eqs. (Ha), (lib), (11c) 
must now be rewritten in terms of moments and slopes. 
This is easily done by substituting Eq. (3b) into Eq. (lib), 
and the result into Eq. (6b) 


where 

A = area of triangle 
v = Poisson's ratio 



M = (T T K„ T) 0 (12a) 


Therefore, the bending portion of the sandwich element 
stiffness matrix is 


A.2 


1 ~f V 
2 


K, = -|-(T r K CJ ,T) (12b) 


Also, it should be remembered that, in Ref. 5, the sheet 
thickness is represented by t, while in this report, h is used 

for cover sheet thickness and t for the total element thick- Carrying out the triple matrix product and using Eqs. (7a) 
ness. and (7b) leads to 




K 


^2 

^2 


^3 ^3 

X“ 

A 23 

+ 





(SYMMETRIC) ] 

- 

\ 2 X:viy-:\ 

A,X 2 3 + 

^23 






* 23 * 3 , 

+ A.t/j.y.,, 

AlX 23 l/l:{ _ 

^, 3 t /23 

*?*+ 





AiXu 

— vX 2A y l:i 

AiX 23 X 3 , + 

t/ 23 l/:u 

A2X13J/31 

AiX“ 3 

4- 

yin 


*12*2:1 

1 d - AjJ/iaJ/a;* 

A 1 X 3 2 J/ 1 2 — 

- tfaif/gs 

* 12*31 d“ AiI/i 2 !/31 

-AiX 1: ,t/i2 - 

VXtat/i;, 

+ ^y?2 1 

1 

1 

"H 

I/23 — vX.vjj ij. 

A1X12X23 d - 

J/iaS/aa 

— AiXi 2 t/l 3 *<Xl 3 j/l 2 

AiX 12*31 

+ 


A2X2 1 1 2 |A,x 2 2 + t/ 2 2 _ 


(12c) 
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For the sandwich beam element, the corresponding flange 
member stiffness matrix can be found from the discussion 
given in Ref. 10. 


D. Stiffness Due to Shear Webs 

It is convenient to begin the discussion by considering 
shear web 1-2 located along the x-axis (Figs. 2 and 7). 
Equation (4.85), Ref. 10, gives the appropriate stiffness 
matrix for this member. Forces X a , X&, etc., are from the 
cover sheets in the case of the sandwich triangular ele- 
ment. Equations (4a) and (4b) apply to the forces shown 
in Fig. 7. 


X« 

Z a 

X, 

Z & 


or 



El— 2 

(SYMMETRIC)” 


t 

t \ 

__ G(A w )i- 2 

-1 

U-2 i 

2 1 

El— 2 

t 

*4-4 

| 1 

i-H 

1 ~ 

1 


1 

L,_ 2 i 1 L 12 J 


X = (K T c)i-2 U 



(13a) 

(13b) 


where {A w ) x - 2 is the cross-sectional area of the shear web 
and the bars over the symbols in Eq. (13b) indicate that 
the member lies along the x-axis. 


The first step is to transform Eq. (13) such that the 
member may lie at an angle ^ to the x-axis. Such trans- 
formations are discussed in detail in Ref. 10. The appro- 
priate transformation equations are 


X = A r X and u — A u 


where 


and 


" A 

/* 

0 

0 

0 

0“ 


A 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

A 

/* 

0 

0 

0 

0 

— fL 

A 

0 

0 

0 

0 

0 

0 

1 _ 


A. = cos 'Jr , ^ = sin * 


(14a) 


(14b) 


(14c) 


It should be noted that, in writing A, Eq. (13a) has been 
rewritten to include Y a and Y a > and v a , tv* This requires 



rows_and corresponding columns of zeros to be included 
in (Km>)i- 2 . This generalization of the equation does not 
alter the basic content of Eq. (13a). Substituting Eq. (13b) 
and the second part of Eq. (14a) into the first part of 
Eq. (14a) gives 

X = A r (K w )j_ 2 A u 

= (K»)i- 2 u (15a) 


Equation (15a) must now be transformed into moments 
and slopes consistent with the terms used for the cover 
sheets. At the same time, Z a of Eq. (13a) should be re- 
placed by Z x (Eq. 4b). Using Eqs. (6a) and (7a), the forces 
of Eq. (15a) may be related to moments as follows (with 
Z x = 2Z a ): 


X a 

Ya 

| Za( 

\x b{ 

Y b 

Z b 


1 MI1 \ 

0 

- 1 

0 

0 

0 

0 

l M Vi 1 

i 

f\ 

V 

A 

u 

A 

V 

0 

0 

) 21 ( 

0 

0 

2 

0 

0 

0 

WrJ ~ 

0 

0 

0 

0 

- 1 

0 

1 ^2 | 

0 

0 

0 

t 

0 

0 

\ Zz 1 

0 

0 

0 

0 

0 

2 

\ / 

— 





— 


(15b) 


By similarly treating displacements (Eq. 3a) 



0 

t/2 

0 

0 

0 

0 

-t/2 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

t/2 

0 

0 

0 

0 

- t/2 

0 

0 

0 

0 

0 

0 

0 

1 
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Now, by substituting Eq. (15c) into (15a) and the result into Eq. (15b) 



(15d) 



Carrying out the triple matrix product of Eq. (15d) gives the shear web stiffness as 


(15e) 




K 

K 

w x 

0* 2 


XV 2 



y'L 



(SYMMETRIC) 




^2 1 1 

* 2 . 







2t/ a , 

2^21 

4 





^ U /1-2 

yh 

— X 2 it/21 

2t/ 2 i 

y 2 , 





*2lt/2l 

X 2 

*21 

2 x 2 i 

^21*/21 

X 2 

*21 





2*21 

-4 


2*21 

4 


(15f) 


The other shear web stiffnesses are simply obtained by 
replacing 1-2 in Eq. (15f) with i — j and then letting i — j 
successively become 1-2, 2-3, and 3-1. The results then 
apply to the triangle arbitrarily oriented in the xt/-plane. 

E. Overall Element Stiffness Matrix 

The total stiffness matrix is now the sum of the con- 
tributions from the cover sheets and three shear webs. 
These are given by Eqs. (12c) and (15f), where the latter 
is written out in full for each of the shear webs. The re- 
sulting stiffness equation is of the form 


/ \ / 

0*i \ 

M Vl 

1 


\ Z ' 

1 

i-i 7 

\m, J 


^2 


= (K(, + (K|«)l-2 “t" (K-U>)2-3 "I" (K«)m) { 

«v 2 ) 

) z 2 ( 


j u; 2 ( 

Im * 3 1 

\ | 

f 1 

I M Wt 


^ 1 

\ Z * 1 

! 

it o 3 J 


(16) 
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Shear flows in Fig. 9b are related to the nodal forces 
of Fig. 9a by putting 

Zi [Q 1-2 h<l-2 Q3-I L3-1] / 2 
Z 2 “ [<7 1-2 f^l-2 "h Q2-3 ^2-3] / 2 

Z 3 “ [<72-3^2-3 ^3-X ^3-l ] /2 

Solving for the shear flows 


on a section at this line be q x . Then, considering forces 
on the portion of the triangle to the left of the cut and 
imposing equilibrium 

qx- 2 * 3 + <73-1 + q x t/ 3 = 0 

Solving for q T and using Eq. (18d) 

q r — 2(Z 2 4- Z% Xa/ L 1 - 2 ) /y * (18e) 


q 1-2 2 1 

<72-3 f^2-3 / 
<73-1 b 3 _! J 



The corresponding transverse shear stress is 

— <7x/t 


(18f) 


The shear forces of Fig. 9a must satisfy equilibrium. 
This is so because these forces have been calculated 
from the known displacements and stiffness equation for 
the element. 

Consequently 


This shear stress is constant for the element. 

Next, consider a section cut parallel to the x-axis. It is 
convenient to select this section at the edge 1-2 as shown 
in Fig. 10b. Then, again imposing equilibrium, and 
letting the shear flow on the section be q y) it is found that 


Z, + Z 2 + z 3 - 0 (18c) 


q» = - qi -2 (18g) 


Substituting Eq. (18c) into Eq. (18b) yields 


and the shear stress by 


q 1-2 = 2Z 3 /L 1 _ 2 

< 72-3 2Z 1 /Zj 2-3 (18d) 

<73-1 — 2Z2/L3-X 

From these shear flows, it is now possible to calculate 

the shear stresses in the element. 

First, consider a section of the solid plate cut normal 
to the x-axis. For convenience, take the perpendicular 
from node 3 to edge 1-2 (Fig. 10a). Let the shear flow 


Tyz - q y /t — 2Z 3 /L,-2# (18h) 

This stress is also constant on any section cut parallel to 
the x-axis. 

3. Alternative procedures . The above stress calcula- 
tions are simple because they follow directly from the 
stiffness method solution. Improved accuracy may be 
achieved by using various refinements on the basic pro- 
cedures. Reference 11 describes some of these refine- 
ments. 


y y 




Fig. 10. Sections in solid element used in determining 
.transverse shear stresses 


III. Examples-Thin Plates and Shells 

A. Flat Plates 

Figure 11a illustrates a flat plate of aspect ratio 2. 
Edges may be either simply supported or clamped. Load- 
ing will be either uniformly distributed over the plate 
area, or concentrated at the midpoint. Figure lib shows 
the elements corresponding to the case n = 2. Because of 
symmetry, only one quadrant of the plate needs to be 
considered in the analysis. For the case n = 4, each 
subrectangle of Fig. lib is divided into four equal 
rectangles, which, in turn, are divided into two triangles, 


1 
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I 


(b) 


Fig. 11. Rectangular plate represented by triangular elements 


as in the case shown in the figure. A similar definition is 
used for other values of n. In all, solutions were obtained 
for n = 1, 2, 4, and 8. 

The displacement at the midpoint may be written as 
Pa 2 

uw = ct c (concentrated load) (19a) 

QQ* 

uw — etd (distributed loading) (19b) 


the convention of Ref. 12, the internal moment (in.-lb/in.) 
acting along the y- axis (Fig. 11a) will be compared 
against the exact result. In Ref. 12, this moment is repre- 
sented by (M x ) x = o, while in the notation of this report 
(Fig. 4a) it is termed (M„) * =0 . Reference 12 gives values 
for where 


and 


(M z ) x =o = Pqa 2 


q is the distributed loading in psi. 


where 


D - Et 3 / 12(l-v 2 ) 


as in plate theory. 

Reference 12 gives the exact solution based on plate 
theory. Table 1 gives the stiffness solutions for the maxi- 
mum deflection. These solutions agree remarkably well 
with the exact solution, even when the idealization of 
elements is quite coarse. 

Bending stresses in the plate are also of interest. Only 
the case of the simply supported plate under a uniformly 
distributed loading will be presented herein. Following 


Table 1. Maximum deflection— rectangular 
plate problems 


n 

Simply supported 

Clamped 

Cte 

a<i 

a c 

a d 

1 

0.0219 

0.0109 

0.0119 

0.0059 

2 

0.0191 

0.0108 

0.0098 

0.0039 

4 

0.0174 

0.0103 

0.0082 

0.0029 

8 

0.0168 

0.0102 

0.0075 

0.0026 

Exact 





solution 

0.0165 

0.0101 

0.0072 

0.0025 


Figure 12 shows /? as taken from Ref. 12. In addition, 
the stiffness solutions, for the idealizations n = 4 and 
n ~ 8, are also shown. To obtain the stiffness solution on 
the y-axis, values for M x were averaged for each pair of 
triangles within a given subrectangle. This average value 
was then taken to hold at the centroid of the subrec- 
tangle. Plotting these results and extrapolating the curves 
to the y-axis then gave the desired results for M x . 



2 y/o 


Fig. 12. Bending moment (M*) along y-axis for 
the uniformly loaded, simply supported 
rectangular plate 
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B. Spherical Cap 

A simple shell is represented by the spherical cap under 
uniform normal pressure. The theoretical solution may be 
found in Ref. 12, pp. 553-554. Figure 13 illustrates the 
problem. 



Fig. 13. Section of spherical cap under uniform 
normal pressure 

The idealizations represented by n = 3 and n = 4 are 
shown in Figs. 14a and 14b, respectively. Other values 
of n are similarly defined. It should be noted that a 
complete quadrant of the cap was used in the stiffness 
analysis. This is basically unnecessary because of the sym- 
metry in this problem; however, the program in use at 
the time the solution was obtained was such that it was 
convenient to use the full quadrant. 

For this problem, the finite element substitution was 
made as illustrated by Fig. 14c. The curved arcs were 
replaced by chords. An interior node was then estab- 
lished by the intersection of straight lines joining oppo- 
site sides of the quadrilateral. Using this interior node, 
each quadrilateral was then represented by four triangles. 
These four triangles were then forced to lie in a single 
plane. This procedure was followed for the subareas 
shown in Figs. 14a and 14b. 


Details of the stiffness method solution in treating this 
problem will not be given here. Within any quadrilateral, 
stresses were averaged for the four triangles and then 
said to apply at the interior node. 

Deflections normal and tangential to the shell are 
shown in Fig. 15. The tangential deflections point radially 
outward from the origin in a plan view. This figure shows 
that if n — 5, the stiffness solution duplicates the deflec- 
tions obtained by exact theory. 



The meridional moments, Mc#>, which again are propor- 
tional to the corresponding internal stresses, are shown 
in Fig. 16. Again, good agreement exists between exact 
and finite element solutions. 



Fig. 14. Idealization of spherical cap into finite triangular elements 
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Fig. 16. Meridional bending moment, spherical cap 


Finally, the so-called hoop stresses, N 0 , are shown in 
Fig. 17. The idealization represented by n — 5 was nec- 
essary to obtain satisfactory agreement with correct re- 
sults in this case. 



Fig. 17. Hoop stress due to bending, spherical cap 


C. Remarks 

The above results illustrate the remarkable ability of 
the flat sandwich element to predict the behavior of thin 
plates and shells. Also, from the structural analyst’s point- 
of-view, it is important to recognize that comparatively 
few elements are needed to achieve results of sufficient 
accuracy to successfully carry out a structural design. 


IV. Orthotropic Panels 

A. Nature of Problem 

Plate and shell structures having orthotropic properties 
are frequently used by the structural designer. A variety 
of design possibilities exist, bringing some confusion into 
the subject. Several basic types of structures will be 
studied herein to show how they may be represented by 
the triangular sandwich element. The discussion will be 
confined to rectangular panels whose sides are parallel 
to the principal directions of orthotropy. 

At the outset, the prime question which arises is the 
following: can the sandwich element discussed in this 
report represent the orthotropic panel to an acceptable 
degree of accuracy? To do so, it must have different 
bending stiffness in two mutually perpendicular direc- 
tions. Before considering this question further, the ortho- 
tropic panels which are to be investigated will be briefly 
defined. 

B. Panels Having Material Orthotropy 

2. Solid plate . This case differs from the conventional 
thin plate because the material properties are orthotropic, 
rather than isotropic. Plywood, in which the grain of 
alternate plys runs at right angles to each other, repre- 
sents a generalization of the basic case. 

2. Sandwich panels . These consist of identical upper 
and lower cover plates separated by a shear core. The 
sandwich triangle will directly represent this case when 
the cover sheets and shear core are isotropic. It is also 
of interest to consider the more general case; i.e., when 
the panel covers are orthotropic. 

C. Panels Having Geometric Orthotropy 

Examples of this type of panel have been suggested, 
tested, and used. Of long standing is the panel made up 
of corrugated sheet held between identical top and bot- 
tom covers. Or, alternatively, the interior structure may 
consist of equally spaced stiffeners such as I-sections, 
channels, or even tubes. Panels of this type are ortho- 
tropic in nature, although the material may be isotropic 
throughout. A generalization occurs when the cover 
sheets are orthotropic. 

D. Orthotropic Stress-Strain Law 

Assume a rectangular sheet for which the two prin- 
cipal directions of material properties are parallel to the 
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x and y edges, respectively. The stress-strain relation for 
sueh a sheet may then be expressed (Ref. 13) as 






v 2 E t 

e 2 

0 


0 " 
0 

2 G 



(20a) 


Symmetry requires that 


v 2 E 1 


Vl £ 2 


(20b) 



Fig. 1 8. Orthotropic plate and idealization into finite 
element triangles 


Equations (20a) and (20b) indicate that four material 
constants are required to define the orthotropic, two- 
dimensional case. However, as for the isotropic material, 
no coupling exists between normal and shearing terms. 

As noted in Ref. 13, the above material properties are 
often usefully replaced by defining 

E — Ely kE — E 2y v — v\y kv — v 2 
G = aE/1 - kv 2 (21) 


Section II-B-2 provided the basis for modeling the 
conventional thin plate with the sandwich element. The 
basic requirement was that the bending stiffness be the 
same for actual plate and sandwich model. A similar 
requirement will now be imposed in modeling the ortho- 
tropic panels by the sandwich element. 

F. Modeling the Solid Orthotropic Plate 

I. Solid single-ply section . 

a. Cover sheets. In the x-prineipal direction put 


The four independent constants are then E, v, k, and A. 
In terms of these constants, Eq. (20a) becomes 


} 

1 F 

“1 

k v 

0" 

\ 

| ~ 1 -kv 1 

kv 

k 

0 

r *v) 

0 

0 

2a 



( 22 ) 


Hence, either Eq. (20a) or Eq. (22) may be used to repre- 
sent the two-dimensional orthotropic sheet. 


E. Modeling an Orthotropic Panel With a Sandwich 
Element 

The orthotropic plate will be taken as shown in Fig. 18. 
A typical choice of triangular elements is also shown. 
This idealization agrees with that used in analyzing flat 
plates in bending (Fig. lib). 

An obvious basis for providing the triangular sandwich 
element with orthotropic bending properties is to use 
orthotropic cover sheets. By assigning the proper values 
to E u E 2 y Gy and Vi (or v 2 ) the stiffness can be altered at 
will in the two perpendicular directions. It would then 
appear that the sandwich element could become capable 
of representing the orthotropic panel. 


Ex It — Ex lx 

x k T a x p x p 

where subscript s stands for the sandwich element and p 
for the true structural panel. As for the isotropic sand- 
wich element (Fig. 5), the moment of inertia I XfS is given by 

Ir H = hd*/2 

And, if t is the thickness of the actual panel 

h p = *712 

By putting 


E Xs = E Xp and d = t 
it is again seen that 


h = t/6 

where h is the thickness of the orthotropic cover sheet on 
the sandwich element. A similar result is obtained in the 
{/-principal direction by choosing 

E v, ~ E v p 
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Consequently, the orthotropic cover sheets on the sand- 
wich element permit the bending stiffness in the principal 
directions to match those of the actual orthotropic plate. 

b. Shear webs. The shear web thicknesses can be de- 
termined by using Eq. (15f) with shearing modulus G 
equal to the value known for the actual plate. 

2. Solid multi-ply section . 

a. Cover sheets . Next, a multi-ply section, similar to 
plywood is considered. The plate consists of lamina, sim- 
ilar to that illustrated in Fig. 19, with a 90-deg change in 
principal directions from one ply to the next. Symmetry 
about a horizontal midplane is assumed. Figure 19 repre- 
sents a unit strip of such a panel section. Any number of 
laminae may be used as long as the symmetry condition 
is maintained. For the actual panel section of n lamina, 
let 

£(£«/*) i =* (^)avm2 

i-1 

where the summation is over the lamina making up the 
panel. This expression can be used to determine (E Xp ) &v . 

Now equate 

Ex s Ixf. (Fxp)av £V 12 

where 

h 8 = hd 2 / 2 

Then, by letting 

d — t and E Xg — 

it is again found that h = t/6. 



Fig. 19. Multi-ply orthotropic plate section 

A similar result is obtained by considering terms in 
the (/-direction. It is then necessary to use 

E>v s = (-Ef/ p )av 

b. Shear webs. The value of G to be used in Eq. (15f) 
is known from the material properties of the orthotropic 
lamina. If the individual laminae have different values of 
G, a weighted average can be used in determining the 
shear web areas of the sandwich element. 

3. Sandwich panel . No difficulty is encountered with 
this type of panel. The sandwich element cover skins can 
directly duplicate the orthotropic covers of the actual 
panel. Also, G of the panel shear core can be used di- 
rectly in Eq. (15f) to determine the sandwich element 
shear webs. The depth, d, of the sandwich element should 
exactly equal that of the actual structural panel. 

4 . Panels having geometric orthotropy. Typical ex- 
amples of panels having geometric orthotropy are illus- 
trated in Fig. 20. Orthotropy exists even when the ma- 
terial is isotropic throughout. 

a. Isotropic cover sheets. Considering the x-direction, 
first calculate (I* ) av . This is the average moment of 
inertia of both cover sheets and interior stiffeners per inch 
(measured along t/-axis). Since (Ix p ) av is known, calculate 
an equivalent solid plate thickness, t , by requiring 

*712 = (h p ) av 


M 



(a) (b) 


Fig. 20. Geometrically orthotropic panels 
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The problem can now be regarded in the same manner 
as the previous solid plate (Subsection F-l). In this way, 
the triangular sandwich element should be modeled on 
the basis of 

h = t/6 

and 

E Xg = E Xp 

where h is again the element cover skin thickness. If E* 
of actual plate covers and stiffeners are not equal, E Xp 
must be replaced by a suitably determined (E Xp ) av . 

A similar calculation can be made in the (/-direction. 
Stiffeners are now ineffective in contributing to I Vp 
(Fig. 20). For equal bending stiffness in the (/-direction 

Ev 8 Iva = Ey p ly p 

For the sandwich triangle, l Vg is the same as I Xg ; i.e., ht 2 / 2. 
Hence, the last equation can be used to find E Vg . This 
completes the determination of the finite element cover 
sheets for this case. 

b. Orthotropic cover sheets . First, calculate l Xp (per 
inch) for plate covers and stiffening material. Then, find 
the equivalent solid plate thickness, t , from 

*712 = l Xp (covers) -f l Xp (stiffeners) 

Then, find (E Xp ) av by requiring that 

(E Xp ) av *712 = E Xp l Xp (covers) + E Xp l Xp (stiffeners) 

This last expression determines (E Xp ) av for the panel. 
The sandwich triangle can now be introduced. For equal 
bending stiffnesses, between sandwich and actual panel 
in the x-direction 

(E x ) g hd 2 /2 = (E Xp ) &v *712 

Letting 

d = * 

E* s = (E Xp ) av 

it is found that 

h - */6 

In the (/-direction, the stiffeners are once again ineffec- 
tive in providing bending stiffness. Hence, I Vp is due 


entirely to the cover sheets. If E Vp is used as the appro- 
priate panel cover sheet modulus, and * and h are retained 
as determined above, equal bending stiffnesses in the 
(/-direction will result if 

Ey g ht 2 /2 E Vp l Vp 

This last equation then determines E Vg for the triangle. 

c. Shear webs. As may be implied from Fig. 20, 
geometrically orthotropic panels will have small shear 
flexibility in comparison to bending. Or, in other terms, 
the strain energy in the triangular element shear webs 
will be small compared to that in the cover sheets. As a 
result, the overall stiffness solution will not be sensitive 
to appreciable errors in determining the shear web thick- 
nesses. The shear webs are necessary but their areas are 
not a sensitive factor insofar as final results are con- 
cerned. This is further discussed in Section VI-C. 

With this in mind the shear web areas can be deter- 
mined by using the equivalent solid plate thickness, *, 
as calculated above, plus G for the actual panel material. 
If covers and stiffeners are of different materials, a reason- 
ably calculated average G should be used in Eq. (15f). 

V. Stiffness Matrix Derivation— Orthotropic 
Cover Sheet 

The discussion in Section IV indicates that a sandwich 
finite element can be used to represent several types of 
orthotropic panels if the element is provided with ortho- 
tropic cover sheets. The key to this development lies in 
finding the stiffness matrix for the orthotropic, two- 
dimensional sheet. Once this stiffness matrix is known, it 
may be used in place of the corresponding isotropic rela- 
tion, Eq. (11c), in deriving the total element stiffness 
matrix. Fortunately, the derivation of the stiffness matrix 
for the orthotropic sheet is no more difficult than for the 
more usual isotropic case. 

Orthotropy is defined by either Eq. (20a) or Eq. (22). 
Axes x and y agree in direction with the principal ma- 
terial directions. Transformation of the above equations 
to arbitrary x', (/' axes is not difficult; however, such 
generality is not necessary when idealizing panels in the 
manner illustrated in Fig. 18. As a result, the derivation 
will be carried out for the case shown in Fig. 21. 

Consider either Eq. (20a) or Eq. (22) to be expressed 
as 

a = $ e (23) 
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V*v | 


i(*3. -K3) 


l£ 


A 





X, u 


1(0,0) 2 (jr 2t 0) 

Fig. 21. Orthotropic triangular cover sheet 

Either equation may be used and the choice made will 
specify The strain energy U in the sheet is then given 
by 

U = -i- jjj a dxdydz = - JJ e T o dx dy 

Substituting Eq. (23) 

u = yJJ S T (.Edxdy (24) 

It is known that the correct choice of displacement 
functions for this case is 

u (x, V) ~ a + bx + cy 
v (x, y) = d + ex + fy 

it follows that 


Writing u (x, y) and v (x, t/) at each node of the triangle, 
and solving the resulting set of six equations for a, b, • • • , 
f yields 

a = Ui 
b = 


U 2 — u x 


x 2 


*3 — *2 

X 3 

ill 

* 2 t /3 

x 2 y 3 

Vi 


V 2 - ©1 


*2 


x 3 - X 2 ^ 

*3 


u 2 H 

?/ 3 


* 2 t/ 3 


t>2 + t>3 

*2*/ 3 ?/3 


Then, since 


du 7 

‘■ = s; = b 

dv _ , 

€v ~ dy~ f 

__ du dv 

y " - Jjj + li 

= c + e 


yM 



— 





— 

1 U 1 \ 

\ 

_ J_ 

0 


0 

0 

0 

l Vi I 

I 

x 2 


X 2 




f 

0 

x 3 ~ X 2 

0 

X 3 

0 

J_ 

) u 2 f 

/ 

x 2 y 3 

x 2 y 3 

ys 

j V 2 / 



1 

x 3 x 2 

1 

* 3 

J_ 

1 

0 

1 Us 1 

/ 

2x 2 y 3 

2x 2 

2x 2 y 3 

2a: 2 

2y 3 

\ V 3 j 


(25a) 


or 


e = 0 u 


(25b) 


where A is the area of the triangle. The form of Eq. (26) 
then provides the following for the orthotropic sheet 
stiffness matrix: 


Equation (25b) can now be substituted into Eq. (24). 
Since the integrand is independent of x and y, the result 
may be written as 


K = hA ® T Z 0 


(27) 


V = ±u T (hA@ T f,®)u 


(26) 


Since © and $ are known, it is now simple to complete 
the calculation for the stiffness matrix of the orthotropic 
sheet. 
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VI. Effect of Triangle Geometry on Element 
Stiffness Matrix 

A. Definition of Terms 

Let 0i be the angle between the two edges of the 
triangle which meet at node i. Then the geometric cases 
of interest are defined by 0* < 90 deg, 0i — 90 deg and 
Oi > 90 deg. In this connection, assume 0i to be the 
largest angle included between two edges of the triangle. 

Stiffness matrix K for the triangular element will be 
symmetric. Consequently, its eigenvalues will always be 
real. Depending on how the stiffness matrix is obtained, 
these values may, however, be all positive, or one or 
more may be zero, or even negative. This last possibility 
signals a defect in the element stiffness matrix. 

By matrix transformation, it is possible to replace K 
with an equivalent diagonal matrix. If this is done, the 
following definitions apply if all elements on the diago- 
nalized form of K are: 

(1) Positive, K is said to be positive definite. 

(2) Positive, with one or more being zero, K is said to 
be positive semi-definite. 

(3) Positive, but with one or more being negative, K is 
said to be indefinite. 

A stiffness matrix which is indefinite is unacceptable. 
The negative element in the diagonalized form indicates 
that, for some applied force, the corresponding displace- 
ment will oppose the direction of the load. This is not 
characteristic of physical systems. Consequently, an ele- 
ment stiffness matrix should always be positive definite, 
or at least positive semi-definite. 

B. Triangle Geometry and Character of Stiffness Matrix 

In a numerical study, Utku and Melosh showed that 
the triangular sandwich element had the following 
property: 

1. positive definite if ^ < 90 deg 

2. positive semi-definite if By = 90 deg 

3. indefinite if > 90 deg 

These results are discussed in Ref. 14. The overall effects 
of this behavior are not known at present. Apparently, 
when (9i > 90 deg, the diagonalized form of K contains 
one negative element; if = 90 deg, one zero term occurs. 
An idealization containing triangles having > 90 deg 
can be somewhat artificially treated so as to avoid numer- 


ical difficulties or poor results. It can be stated that, hav- 
ing recognized the shortcoming, steps should be taken to 
avoid it by somehow circumventing the development of 
an indefinite form for the element K matrix. Utku and 
Melosh suggest several ways of doing this (Ref. 14). The 
scheme used in the COSMOS and SAMIS computer pro- 
grams is to use absolute values on certain minor terms 
which are responsible for the generation of an indefinite 
stiffness matrix. Although indefensible on theoretical 
grounds, experience shows that this artifice works. For 
example, the idealization shown in Fig. 14 contains two 
triangles, in each quadrilateral, for which 0* > 90 deg. 
Yet, by using the scheme previously mentioned, excellent 
data have been obtained. As shown in Section III-B, 
extensive experience with the triangular sandwich 
element supports this conclusion. 

The reason for the generation of the indefinite K matrix 
will be discussed in this section of the report. For certain 
classes of triangles, having 0i > 90 deg, an alternative 
suggestion will be made to avoid the indefinite stiffness 
matrix. The suggestion will be based on a logical special- 
ization of the theory which has been used in deriving the 
element K matrix. 

A useful idealization for rectangular planforms is illus- 
trated in Fig. 22. Each such subrectangle is divided into 
four triangles. This automatically introduces two tri- 
angles for which > 90 deg. 

The idealization illustrated in Fig. 22 is essentially the 
same as that used in analyzing the spherical cap of Fig. 14. 
Whether the triangles are plane stress elements used for 
analyzing a two-dimensional continuum, or sandwich 
elements for representing plate or shell bending, stresses 
are always obtained with superior accuracy if the ideali- 
zation shown in Fig. 22 is adopted. The stress at the 
centroid of each subrectangle is taken as the average of 
the stress values for each of the component triangles. 
Because of this, triangles of the type shown in Fig. 22, 
and with 0i > 90 deg, become of particular importance. 
They will be given special attention in this report. 


Fig. 22. Idealization of rectangular forms into four 
triangles 
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Fig. 23. Triangle with 0 3 > 90 deg 


C. Indefinite Form of Triangular Element Stiffness Matrix 

As previously mentioned, Utku and Melosh discovered 
by numerical calculations that K for the sandwich ele- 
ment becomes indefinite when 0* > 90 deg. To examine 
this on the basis of the derivation given in this report, 
the single element k 9y ^ * will be written for a triangle 
taken as shown in Fig. 23. 

From Eqs. (12c) and (15f), where the latter must be 
written for all shear webs, it can be shown that 


- -^r( A * + «&) + -j- (if) i2 *3, ] 


where 


<f> = %(1 -v 2 )A 
Ai = (l-r)/2 
Ate = shear web area 

Xij = Xi — Xj . 


The shear web areas are obtained from Eq. (10) as 

( t) m ~ ~ 8 A ( '~ L '- 2 + L 3- 3 + L '-> ) 

(~ ) 1 _ 3 ~ 8 A ( L ’~ 2 + L '- 3 _ L i - 3 ^ 

Applying the Law of Cosines to the triangle in Fig. 23 
Lj_ 2 Lj 3 + L 2 3 2Li-3 1*2-3 COS 0 3 

Lj_ 3 = Lj_ 2 '+ L 2 3 — 2L x _2 L 2 -3 cos 0 2 
Substituting into the expression for <> y ^ yields 

k ° y ' '*1 ~ 16(1 + v)A £ 3(1 - 7 ) ( y " + T - * 23 ) 

H 2 {(1*1-8 c°s S 3 ) x 2 , + (1*1-2 cos #2) x 2 j } J 


Now, a special case is taken in which 0i = 0 2 = 30 deg 
(Fig. 23). Substituting this geometry then leads to 


k = Et j 2 /. 2 \ 

'*1 9 vi 64(1 + V )A ■ « y 4 13 ) 


Since t is the thickness of the sandwich element, kO Vl 
will be negative. Hence, an applied moment M Vl will 
cause a negative e yl to occur. Due to symmetry, a similar 
condition will exist at node 2. This is physically unrea- 
sonable. It is this unexpected result which underlies the 
indefinite form of K, as reported in Ref. 14. 

It is now of interest to check the limiting case; i.e., the 
triangle for which 0i = 0 2 — 45 deg, 0 3 = 90 deg. If this 
is done 


V 16(1 + v)A 

x hi± r * ( 1 + + -^1 

2 |_3(1 v) y yj2 ) J 


which is always greater than zero. For this same case 


k — — L 2 

4A 2- 


Hence, the indefinite form of K is clearly dependent 
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Further Insight into the problem is obtained by exam- 
ining the shear web areas as given by Eq. (10). From 
Eq. (10) 


(t) 1 . 2 = ^(- L - +L - +L -) 

( t ) 2 _3 = 8 A ( L '~ 2 ~ L *- 3 + 

(t) 3 _, = 8 A ^ 1 ' 2 + L *- 3 ~ 


Substituting from the previous expression for L\ 2 



(2Li_3 Li 2 -3 cos 0 3 ) 


When 0 3 > 90 deg, shear web area (A w )!_ 2 becomes nega- 
tive. The other two shear web areas remain positive. 
This result is general in that it is always the shear web 
area on the edge which is opposite to the obtuse angle 
which becomes negative. If 63 — 90 deg, (A tt> )i_ 2 becomes 
equal to zero. For -0 3 < 90 deg, all shear web areas are 
positive. Even when (A w )!_ 2 becomes negative, the total 
strain energy in shear, Eq. (9), remains positive. For 
example, in Eq. (9) 

(Ato/L)i_ 2 + (A mj /L) 1 _ 3 


or 


There is no known basis for completely overcoming 
this difficulty. In Section D below, an approximate cor- 
rection is offered. This correction takes the shear energy 
of Eq. (9) into account. 

D. Suggested Basis for Avoiding the indefinite Form of K 

The numerical work of Utku and Melosh (Ref. 14) 
and the above discussion indicates that the indefinite 
form of K is caused by the negative shear web area. 
Consequently, a basis for avoiding negative shear web 
areas will be proposed. The discussion will be centered 
on isosceles triangles of the type illustrated in Figs. 22 
and 23 with = 0 2 . 

When 0 3 < 90 deg, no problem exists; all shear web 
areas are positive. When 0 3 = 90 deg, (Au?)i -2 becomes 
equal to zero. Thereafter, as 0 3 increases beyond 90 deg, 
(A w )i - 2 becomes negative. It is this latter case which is of 
interest. 

Equations (8) and (9) represent the transverse shear 
energy in the solid and sandwich elements, respectively. 
Equating these energies leads to the shear web areas as 
given by Eq. (10). The negative area then follows as soon 
as 0 3 > 90 deg. No general basis for avoiding this situa- 
tion seems to exist. 

Consider now the case for which 0 3 = 90 deg. Shear 
web area (A ^-2 is then zero. Rather than permit this to 
become negative as 0 3 increases beyond 90 deg, impose 
the restriction that 


(A w /L)i_ 2 "I - (Ato/Lt) 2 -; 


(Atp)i- 2 0 


will always be positive. This is verified by substituting 
the above expressions for ( A w /L)i- it 

It is the negative shear web area which is responsible 
for the indefinite K matrix. The scheme previously men- 
tioned avoids this difficulty by using absolute values for 
shear web areas. If this is done, K never becomes indefi- 
nite. However, the energy balance leading to Eq. (10) is 
then upset. This still leads to an acceptable stiffness 
matrix because this energy is only a small part of the 
total energy in the strained element. Consequently, it can 
be “tampered with” without seriously impairing the 
ability of the element stiffness matrix for representing 
bending behavior. However, if a sandwich-type panel 
having considerable shear flexibility is to be analyzed, it 
is not certain that the absolute value scheme for shear 
web areas would lead to satisfactory results. 


when 


0 3 > 90 deg (28a) 

Further, impose the restriction that 

(Atf>)l-3 = (A«j) 2 -3 — A w 


when 


0 3 > 90 deg (28b) 

Finally, replace the general deformation modes 
(wx^Wt^iVs) in Eqs. (8) and (9) with the average 
condition that 


U>! = u> 2 = ttf 3 = tu 


(28c) 
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Under these restrictive assumptions, the energy equality 
between Eqs. (8) and (9) leads to 


A w — 


t 

8 A 


/ LU + L»_ 3 + Lt A 

\ E t -3 + L 2 -3 / 


^ 1-3 


L 


2-3 


(29) 


This expression can then be used for the isosceles triangle 
with 0 3 > 90 deg. 

In applying Eq. (29) to panels having considerable 
shear flexibility, it may be worthwhile to introduce a 
numerical factor a such that 


A* — a A w (3^) 

The parameter a can then be assigned a numerical value 
such that the shear web energy under the above restric- 
tions represents the best possible average for a range of 
values of 0 3 > 90 deg. For conventional thin plates and 
shells, this is believed unnecessary due to the relative 
unimportance of the shear web terms. 


Numerical data have not yet been secured with Eq. (29). 

VII. Large Deflection and Stability Analysis 

As mentioned in Section I, the triangular sandwich 
element can be used in carrying out large deflection and 
stability calculations on plates. It may also be applicable 
in this same sense to shells, although actual calculations 
have not yet been carried out. 

The reason for this useful characteristic of the sand- 
wich element is that it remains a flat surface while 
representing bending behavior. As a result, it remains 
compatible with the initial stress stiffness matrix devel- 
oped by the author for a triangular element (Ref. 9). The 
ability of these two stiffness matrices for the triangle 
(sandwich element matrix of this report and the initial 
stress matrix of Ref. 9) to correctly predict plate stability 
and large deflection behavior was shown in Ref. 15. This 
opens the door for a finite element solution for many 
types of plates, in which stability or large deflections, as 
defined by the von Karman equations is of interest. 
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